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Abstract 

We study the behavior of light rays in perturbed Robertson- Walker cosmologies, cal- 
culating the redshift between an observer and the surface of last scattering to second order 
in the metric perturbation. At first order we recover the classic results of Sachs and Wolfe, 
and at second order we delineate the various new effects which appear; there is no a priori 
guarantee that these effects are significantly smaller than those at first order, since there 
are large length scales in the problem which could lead to sizable prefactors. We find that 
second order terms of potential observational interest may be interpreted as transverse and 
longitudinal lensing by foreground density perturbations, and a correction to the integrated 
Sachs- Wolfe effect. 
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I. Introduction 

In the last several years, observations of temperature anisotropies in the cosmic mi- 
crowave background (CMB) [1] have spurred increasingly sophisticated investigation of the 
anisotropy predicted by theoretical models [2-5] . Important contributions to the anisotropy 
come from gravitational perturbations, temperature and pressure fluctuations at the sur- 
face of last scattering, and ionization effects in the later universe. 

The earliest of these effects to be studied, and the most important on large scales, 
are those due to gravitational perturbations. These were systematically investigated by 
Sachs and Wolfe [6], who derived the basic formulae relating perturbations in the metric 
to anisotropy in the temperature of the CMB. Their results revealed two basic sources 
of anisotropy: potential fluctuations at the surface of last scattering, and time variation 
of the potential along the path of the photon. Later investigations focused on individual 
effects in speciflc models [7-15]. 

Even though perturbations in the energy density 5p/ p grow to be greater than unity on 
sufficiently small scales, the resulting metric perturbations may almost always be taken to 
be small [16]. It therefore makes sense to calculate the behavior of photons to first order 
in this pertubation, as Sachs and Wolfe did. Nevertheless, there is no way of knowing 
ahead of time that second-order terms in an expansion in the metric perturbation will be 
negligible compared to the first order terms, since there is ample opportunity for effects 
to accumulate as photons travel to the observer from the surface of last scattering; in 
other words, the coefficients of the second-order terms may be numerically large. (As an 
example of a related effect, the time delay formula in standard gravitational lens systems 
contains important contributions from both the first order Shapiro and the second order 
geometric effects.) It is therefore worthwhile investigating the redshift induced by effects 
which arc formally second order in the metric perturbation to see if they may nevertheless 
be observationally important. In this paper we calculate these second-order effects and 
interpret the results in terms of specific physical processes. 

It is necessary to be careful about what we mean by "second order" in the context 
of gravitational perturbation theory. We imagine that we are given a metric throughout 
spacetime of the form 

= 9^^y + V ' (1-1) 
where g^^u describes a background Robertson- Walker spacetime and h^i, is a small per- 
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turbation. We will not be computing this perturbation to second order in pertubations 
of the energy-momentum tensor, but simply calculating photon trajectories to second or- 
der in hf^u and its derivatives. Therefore, if h/j^i, is computed from standard first-order 
metric perturbation theory and substituted into our expressions, the results will not repre- 
sent a complete calculation of effects which are second order in the matter perturbations. 
(In Sec. IV we will examine explicitly the case of first-order scalar perturbations, but it 
is straightforward to generalize the results.) Nevertheless the expressions we obtain will 
constitute a subset of all the possible contributions, and if any of them turn out to be 
comparable in magnitude to terms which are formally first order, it is appropriate to take 
them into account. Moreover, the substitution h^i, i— g^^J + gj^J + . . . into our formulae 
below would immediately yield an expansion for the full second order anisotropy. 

It is also important to note that we will only be dealing with gravitational perturba- 
tions. We will imagine that there is a hypersurface of last scattering fixed at some definite 
time, on which there can exist intrinsic perturbations which may be calculated indepen- 
dently; we then compute the additional perturbations due to the metric fluctuations along 
the geodesies followed by the photons. Non-gravitational second-order perturbations were 
treated by Vishniac [17], Dodelson and Jubas [18], and Hu, Scott and Silk [19]. The lat- 
ter authors also examined higher order gravitational effects by expanding the Boltzmann 
equation to second order, but did not construct explicit solutions. As a final caveat, we 
treat the order-by-order expansion in powers of the metric perturbation and its derivatives 
in a formal sense; thus, a phenomenon such as the integrated Sachs- Wolfe (or Rees-Sciama) 
effect we consider to be first order (since it involves terms linear in derivatives of /i^jy), 
even though it is sometimes thought of as second order since it can be numerically small 
(and vanishes to first order in some specific models). 

Our calculation proceeds as follows. In Sec. II we set up the problem and express the 
redshift experienced by a photon in terms of its corresponding background path x^'^^^{X) 
and its first and second order perturbations, x^^^^{X) and x^'^^^{X). In Sec. Ill we discuss a 
general formalism for constructing these pertubations in terms of the metric variables; this 
is an extension of the methods of Pyne and Birkinshaw [20] to arbitrary order. In Sec. IV 
we specialize to the case of scalar pertubations, and examine the resulting formula for the 
temperature anisotropy. Although a quantitative understanding of the magnitude of each 
term would require detailed knowledge of the evolution of the metric perturbations (which 
we do not attempt in this paper), it is possible to discuss informally which contributions 
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might be observable in realistic models of structure formation. 

II. Perturbation Expansion 

We are interested in the pattern of temperature fluctuations AT/T on the sky as seen 
by an observer in a perturbed Robertson- Walker spacetime. We write our background 
metric in conformal coordinates = {rj, x, y, z) as 

ds2(0) = g<;^Jdx^dx'' 

(2.1) 

= a'^{r])[-dTf + (dx^ + dy"^ + dz^)] . 

Here 7 = l + nr'^/A, where n is the spatial curvature parameter (+1, —1 and for positively 
curved, negatively curved and flat cases, respectively), 0(77) is the scale factor, and = 
x^ + y"^ + z^. In this section we consider an arbitrary metric perturbation h^i,. It will be 
convenient to separate out the dependence on the scale factor by working in the conformal 
background metric gj^J = a~'^g^^u, with the conformally-transformed perturbation h^^j ~ 
(i~'^h^iy, so that the actual, physical spacetime metric is given by g^i, = gl^J + h^u- The 
wavevector of a light ray in the physical metric is related to the wavevector in the 
conformally transformed metric by k^^ = a^k^. (Our conventions are those of Ref. [20].) 

Within such a spacetime we consider a photon path a;'^(A), where A is an affine pa- 
rameter. (See Fig. 1.) This path connects an observer at a point O with coordinates 
Xq = (770,0,0,0) to the hypersurface of emission, which we deflne to be the spacelike 
hypersurface of constant conformal time r] = rjs- The "surface" of emission is then the 
intersection of the past light cone of the observer with this hypersurface. We assume that 
at conformal time ryg every point with spatial coordinates emits thermal radiation with 
a temperature T£-(p, d), as measured by a comoving observer, which depends both on po- 
sition and on direction as characterized by a three- vector d, normalized to imity in the 
background metric, gjj^J , restricted to the hypersurface. (This hypersurface need not be 
the actual time of last scattering, but need only represent a hypersurface on which the 
radiation fleld is understood.) The photon path itself is specifled by a three-vector e in the 
hypersurface of constant conformal time containing O normalized to unity in gjj^J . We can 
think of e as the direction on the sky toward which a comoving observer at O is pointing 
an antenna; for observers which are not comoving e and the observer's direction vector are 
related by a Lorentz transformation. The initial condition e determines the point p and 
direction vector d at which the ray intersects the hypersurface of emission. 
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To an observer with four- velocity U^^ (normalized to U^^g^i^W^ = —1), a photon with 
wavevector k'^ = dx'^/dX, with A an aflfine parameter in the conformal metric g^^, has a 
relative frequency given by 

to = -a-^g^^U^'k" . (2.2) 

(We refer to this as the "relative" frequency, since we are free to scale the affine parameter 
A to set the normalization of oj. The ratio of relative frequencies at two points along the 
path is invariant under such a reparameterization.) For a blackbody spectrum, the CMB 
temperature observed at O is related to the temperature at emission by 

To(xo,e) = ^T£(p,d) . (2.3) 
Us 

We are therefore interested in computing, given the initial data x^,e, wo, the quantities 
p, d, and uj£. These depend on the photon path and associated wavevector, which we may 
express as series expansions in the perturbation h^i, and its derivatives: 

x^(X) = x(°)^(A) + x(i)^(A) + x(2)^*(A) + . . . 

(2.4) 

k^'{X) = /c(°)'^(A) + /c(^)'^(A) + A;(2)'^(A) + . . . 

The situation is thus as portrayed in Fig. 1. Note that x^^^^{X) has the interpretation of a 
path through spacctime, while the x^'^^^{X) arc thought of as deviation vectors at each A. 
In this section wc will calculate the observed temperature in terms of these quantities (plus 
the intrinsic temperature fluctuations on the surface of emission) , while in the next section 
we will explicitly calculate the path and wavevector in terms of the metric perturbation. 

We have already specified Tg as the temperature measured by a comoving observer. It 
will also be convenient to take our observer at O comoving. This requirement is physically 
acceptable, since any motion of the observer leads to a dipole anisotropy which may be 
easily subtracted. It is sometimes useful to imagine a family of comoving observers with 
four-velocity defined over all of spacetime. The normalization condition U^gf^i^U^ — —1 
then leads to 

[/(O)/^ = 0-1(1,0,0,0) 

U^'^^^ = a-^ Q/ioo,0,0,o) (2.5) 
f/(')'^ = a-lQ(/^oo)^0,0,o) . 



(2.6) 



We can also explicitly construct the geodesies of the background metric, x^^^^^{X). We 
consider null rays which intersect the observer at the spatial origin of co-ordinates, and we 
choose the affine parameter such that 

= 1 , 

g^fk^'H^'^^ = 1 . 

A two parameter family of such rays which satisfy these conditions is given by [21] 

^(o)m = (A,reO 

where the e* are components of e, and 

Xo — A 



(2.7) 



r(A) = 2tanK 
7(A) = sec^ 



2 

Xo — X 



(2.8) 



where Xq is the affine parameter at the observer. The subscript k on a trigonometric 
function denotes a set of three functions: for k = 1 the trigonometric function itself, for 
K — —1 the corresponding hyperbolic function, and for k = the first term in the series 
expansion of the function. (Thus, sIuq 6* — 9, cosq 6* = 1.) Finally, we can place boundary 
conditions on the higher-order quantities x^^^^^, x^'^^^^, k^^^^ and k^"^^^ at the origin. For 
convenience we will set 

xW'^(Ao)=x(2V(Aa) = 

k^^^' (Xo) = k^^^' (Xo) = . 
Then the condition that the wavevector be null at the observer implies that 



(2.9) 



o 



k^^^\Xo) 



[{hoof + hoohoik^'^' + ^/ioo%A;(o)^A;(o)^' + ^ihoik^^^'f 



o 



(2.10) 



o 



The temperature at emission can be written as a uniform background plus a small 
perturbation, expressed as 



Tf(p,d) = [l + T(p,d)]T^' 
6 



(0) 



(2.11) 



The function r will be treated as first order {i.e. of the same order as /i^uz^), and will be 
unspecified in this paper since our interest is in the gravitational effects on photons in 
the time since emission. The point at which the geodesic intersects the surface r] = r]£ 
can be written as p = p^°) + p^-*^) + . . . (Note the distinction between a;*, the spacelike 
components of the separation vector, and p\ the separation of the intersection points of 
the path at different orders with the constant-time hypersurface.) Expanding los and d as 
well, eq. (2.3) to second order becomes 

To = + r(p(0) + p(^), d(0) + d(^))]Tf . (2.12) 

With the conventions chosen in the previous paragraph, lVq^ — a (rjo) ^ and lv^^ = 
a{'qs) ^- The quantity of interest to us is the fractional deviation in the observed tem- 
perature with respect to the expected temperature in the unperturbed spacetime, and we 
denote this deviation by 5T. Expanding r in a Taylor series, we obtain 



5T = 



JP \ To 



, ,(0) / rp{0) 



l+(4"-4'>+r) 

+ (4^' - 4' + - 4)4" + 4"r - 4V + <;">'^) 

■(2.13) 

where the are the components of d, r and its first partial derivatives are evaluated at 
(p^°\d^°)), and we have put uj^°'^ = u!^°'^ / u'^^K We note that our freedom to choose T^°^ 
may be used to render ST observable, e.g. by setting tJ°^ = a (770) a ivs)~^ (To) where 
the angle brackets denote an average over the observer's sky. 

Expanding the metric perturbation and photon wavevector around their values on the 
background path, we obtain 

cl;(°) = 1 

2 ox^ ax^ dX dX 
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In this expression AA is the difference in aflfine parameter between the point where the 
zeroth and first order geodesies intersect the hypersurface 77 = constant; to this order 
AA = -a;(^)°. It is also straightforward to show that p^^^* = x^^^' - A;(°)*a;(^)°, and that 
is given by 



(2.15) 



where the norm is defined by the spacehke part of the background metric. Putting it all 
together we obtain 



(0) 



6T 



5T = 



+ 



o 



o 



+ 



+ 



l;,..J.(0)il.(0)j 



O 



o Z Z 



Z dX 



dd^ 



(2.16) 

Here, the notation £ means that the quantities referred to should be evaluated at the point 
ivs^P^^^) ^iid direction d^^^. 

To complete the above formulae, we have to solve for the perturbed geodesies at first 
and second order in terms of h^^,^. In the next section we carry this out for arbitrary metric 
perturbations, and in the following section we specialize to scalar perturbations. 



III. Second-Order Geodesies 

In order to calculate the approximate geodesies of ^^j^ = g'^J + h^i, order by order 
we employ the perturbative geodesic expansion introduced in Pyne and Birkinshaw [20]. 
Because those authors worked only to first order it is necessary slightly to extend the 
equations to address the higher order questions we are concerned with here. In this section 
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we describe the needed extension, which writes a general solution for the approximate path 
at any order without restriction on the perturbed spacetime under consideration. In the 
following section we specialize this general solution to gain the null geodesies to second 
order of perturbed FRW spacetimes in the longitudinal gauge. 

We begin with the geodesic equation in the metric g/^i, = gj^J + h/^,^, 

^+r^^f,k^k^ = o, (3.1) 

which holds along some path x'^(A). We seek to approximate that path to any given order 
by solving for the x'^°'^^{X) in (2.4). To this end we substitute (2.4) and the equation 

a/3 + r^^^'^a^ + r^^^'^Q./j + ... (3.2) 

into (3.1) and simultaneously Taylor expand each of the F'^'^^^q,^ at x^^{X) about their value 
at x^^^^{X). In (3.2), F^^^'^q./j is that part of Ti^af3 which is of a-th order in either hn^, its 
first partial derivatives, or their products. The resulting equation, equivalent to (3.1) but 
holding along the path a;(°)'^(A), is written 



^a=0 



(3.3) 

At zeroth order we find that a;*'^-''^(A) is an affinely parametrized geodesic in the metric 
gj^J. At every order above zeroth equation (3.3) may be rearranged into the form of a 
forced Jacobi equation for the a-th order separation vector, x^°'^^^{X): 

+ 2F(°)^«^A;(o)«A;(«)^ + a,F(o)^«/3A:^°^"A;^°^^a;(«)'^ = /(«)'^ . (3.4) 

Importantly, the highest order x*-^-*^ or appearing in /(")^ is of (a — l)-th order. For 
instance, the forcing vectors at first and second order are given by 

= _r(l)^^^;,(0)a;^(0)/3 

- ..F(^)^.,.(^)^.(0)«.(0). - ^....F(0)..,.(.^.(.^,(0)«.(0). ^'-'^ 
-F(2)m^^/,(0)«/,(0)^ , 



Pyne and Birkinshaw [20] showed how a general solution to equation (3.4) could be 
written down in terms of the parallel and Jacobi propagators of the background metric 
g^^v ■ These propagators are matrix- valued functions of a pair of points connected by a 
geodesic, and are defined by path-ordered exponentials along the appropriate geodesic. 
The parallel propagator P(Ai, \2Yu is a 4 x 4 matrix given by 

P{X2Ai) = Ve-KY>(^\jjA{X)d>^ , (3.6) 

where V denotes the path ordering symbol and ^ is a 4 x 4 matrix defined hy A^^,^ = 
2/j(0)ap(0)M^^_ The parallel propagator lives up to its name, in that P (Ai, A2)'' ^.v^ (A2) 
is the vector obtained by parallel propagating v^^ from A2 to Ai along the geodesic. The 
Jacobi propagator is an 8 x 8 matrix given by 

t/(A2,Ai) = Pexp(^|^^ (p(Al,A)7^(A)P(A,A0 j) ' ^^'^^ 

where ^Z{X)^'a denotes the 4 x 4 matrix R^°^^'„p^k^^^''k^^^'' evaluated at x^^\X), and and 
1 denote the 4x4 zero and identity matrices, respectively. The Jacobi propagator serves as 
a Green's function for the Jacobi equation in the background spacetime. More information 
about these objects can be found in [20, 22]. 

The solution for a;(")'^(A) and /c(")'^(A) at some affine parameter A2 can now be obtained 
from their values at some fixed affine parameter Ai via 

[6.1 



L ^^^^'^Kj'(Ai,AU(«)(A))^^ 



the integral being taken over the zeroth order geodesic, x^^^^{X). The program for recursive 
calculation of the x^°'^^(X) is now established: having obtained x^°'~^^'^(X) we can solve for 
j(a)^ and thus obtain x^°'^^(X) from (3.8). The recursion starts by solving the geodesic 
equation of the background for some a;^°^'^(A) and then calculating its associated parallel 
and Jacobi propagators. 

The parallel and Jacobi propagators for the radial, null geodesies of g^^J , the confor- 
mally transformed Robertson- Walker metric, were obtained in [23] and are written 

P{X2,Xir.= (!.i 7(A°/).») (3.9) 
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and 

U{X2,Xi) 



cos^(A2-Ai)J sin^ (A2 - Ai) J 
-Ksin^ (A2 - Ai) J cos^(A2-Ai)J 



+ 



respectively. In (3.10) J is a 4 x 4 matrix given by 

0, 



0* - e'cj 



(1-J) (A2-Ai)(l-J) 
(1-J) 

(3.10) 



(3.11) 



(3.12) 



Given a background geodesic specified by the direction cosines e*, any three- vector v*(A) 
may be decomposed into the sum of a longitudinal part (A) pointing along the geodesic 
and a transverse part v]_{X) perpendicular to the geodesic (in the spacelike hypersurface 
of the background), where 

^(A) = eV,V(A), 

Thus, the matrix J serves to project a four-vector into the plane transverse to the photon 
direction in the comoving spatial hypersurfaces. 

The relatively simple form of these propagators for the case of a Robertson- Walker 
metric allows us to obtain the perturbed geodesic and wavevector from (3.8) immediately. 
Imposing the boundary conditions (2.9) at the observer, we obtain 

X 



(«)0(A) = (A - Ao)/c(«)°(Ao) + / (A - A')/^")°(A')rfA' , 
x\^^\X) = 7(A) /' (A - A')7-^(A')/|{"^A')rfA' , 
x'f^'iX) = 7(A) sin^A - X')^-\X')f[^^\X')dX' , 

/•A 

A;(«)0(A) = k^^^^Xo) + / /^")°(A')ciA' , 
^i"^A) = 7(A) / 



(3.13) 



/c?^A) = 7(A) 



Ao 
A 



nr{X) 



(A - A') 



^-\X')&^\X')dX' , 



cos,(A-A')-^sin«(A-A') 



^-\X')f^^^\X')dX' 



These expressions are valid for any metric perturbation; the specific perturbation is en- 
coded in the vectors f^°■^^^. Note that we have written the integrals as proceeding backwards 
along the path from the observer to the point A on the background path. 
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IV. Scalar Perturbations in the Longitudinal Gauge 

In this section we carry out the program described above to second order for the 
metric perturbation /i^j^ given by 

h^^dx'^dx'' = -2(i)drl^ - 2'ii;-f-^ {dx^ + dy^ + dz^) (4.1) 

describing scalar perturbations in the longitudinal gauge [24] . We note that in this gauge, 
(j) and ij) coincide with the gauge invariant metric variables of [24], $ and This will 
allow us to obtain gauge invariant expressions for the observables of interest by replacing 
(j) with $ and ij) with in our final formulae. Of course, only the first order expressions 
are rendered gauge invariant because $ and ^ are themselves gauge invariant only to first 
order. 

To compute the perturbation vectors /f^")/^ we need to calculate the Christoffel symbols 
to various orders. These are given by 

r(°)°.. = 
r(«)^oa = 

^''''^k = -l-^{5^kxf^5,,xf-5,,xf^) 

r(i)o..-_MA.. 

= -doijSij (4.2) 

r(2)Oo^ = -20a«0 
(2)0 _ X 

i 13-^2 



Calculation of the first order vector f^^^^^ proceeds straightforwardly, using the nor- 
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malization g^fk^^^^k^^^^ = SijC^e^ = 1. We find that 

/«° = ao(0 + V')-2^, 

According to (2.16), the only second-order quantity (as distinguished from products of 

— ^(2) 

first-order quantitites) which enters the formula for ST is the timelike component of 
the wavevector, k^"^^^. We therefore do not need to calculate the entire second-order force 
vector, but only the timelike component. In doing so we make use of the decomposition of 
the directional derivative of a scalar along the path into partial derivatives, 

^ = ^00 + k^'^^d.cp (4.4) 

(since k^^^^ = 1). Another relatively straightforward calculation yields 

/(2)0 = _2 A (^k^DO^ + a;(i)-a.</.) + 2A;(^)°ao(</. - V) 

+ a;(i)"9^5o(0 + ^) + 2{(j) + ^)do^ . 

Substituting (4.3) into (3.13), we obtain the first order perturbed geodesic: 

a;(i)o (A) = (A - Ao) {4> - V')o + / [-2</> + (A - A') do{tl^ + </>)] dX' , 

J Ao 

x^^\X) = (A - Ao)(<^ - V')oA;(o)^(A) + /c(°)^(A) / [{^P -</>) + (A - A') do{tl^ + </>)] dX' , 

Ao 

/•A 

x^l^\X) = 7(A) / sin« (A - A') ^{X') [e'e^ - 5'^] djicj) + i;) dX' . 

Ao 

(4.6) 

The explicit construction for the wavevector perturbation, k^^\X)^ may be obtained 
either by differentiation of (4.6) above or from (3.13) directly. In either case 

(A) = ((/>- V')o - 2(/>A - /isw(A) , 

^ (A) + /,(o)^ (A) [(0 - V;)o - (0 - V)a - /isw(A)] , 



(4.5) 



k^l^' (A) = 7(A) / 

J Ao 



2 

^ ^ Kr{X) . 



cos«(A - A') — sin^(A - A') 



7(A') [e'e' -8'^] djicf) + ilj)dX' , 

(4.7) 
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where the integral 

/isw(A) = - [ do{i; + (P)dX' (4.8) 

represents the conventional integrated Sachs- Wolfe effect. 

As noted above, to compute the second-order effect we need only the time component 
of the second-order wavevector; this is given by 

ki^)o^X) = -1(0 + - 2{x^'^^d^4> + k^^^»x + /2(A) , (4.9) 



where the integral I2 is defined as 

-A 



/•A 

(A)= / 2A;WOao(<^_^) + 2((/, + V')ao^ + xW''apao((/) + V) (iA' . (4.10) 



Having obtained the perturbed geodesies and wavevectors in the longitudinal gauge, 
it remains only to substitute into (2.16) to obtain final expressions for the temperature 
anisotropy. At first order we recover the conventional Sachs- Wolfe result, 

5T^^^ = (0 + r)s -(l>o + IiswiXs) , (4.11) 

where once again the notation £ means that quantities are evaluated at the position and 
direction of the intersection of the background geodesic with the surface of emission. The 
second order anisotropy, the main result of this paper, is given by 

. ,(2) 3 1 

=^(l>o - (l>o(f>£ - (f>or£ - -(I)} + (psTg 

- {2(f)o - ipo - (l>£ - Ts - /isw(Af )) /isw(Af ) - hi^s) (4.12) 

+ {x^l^' + /TDfc(°)^) _ d,{cp + t)s + x(i)°ao(</> + V')5- + , 

where the integral 

/td(A)= I {(j) + i;)dX' (4.13) 

is the Shapiro time delay along the path. 

An accurate appraisal of the magnitudes of the various terms contained in (4.12) 
would require knowledge of the initial conditions and evolution of the perturbations 
and i/j, including nonlinear effects. This information is model-dependent, and we will not 
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attempt such a task here. It is nevertheless possible to remark on the possible importance 
of the different effects to observations, based simply on the form in which they appear. 

The quantities 0o, i/jo, and /isw(A^) are all small (< 10~^) in conventional 

models of structure formation. Therefore the terms in (4.12) which are written as products 
of these numbers are even smaller, and should not contribute to the anisotropy at an 
observable level. Similarly the term S^^'^{dT/dd^) will typically be the product of two 
small quantities, and may be neglected. Therefore the potentially interesting terms are 
those involving the separation vector a;^^^'^(Ag) (which is not necessarily small) and the 
integrals Itd{^£) and /2(Af ). 

The term x^^^^di{(j) + t)^ is due to the transverse deflection of the photons by sources 
between us and the surface of emission; this effect has been studied previously in investi- 
gation of the impact of gravitational lenses on CMB anisotropy [25-36]. (The processing 
of CMB anisotropy by lensing is second order since both the lens angle and the initial fluc- 
tuations being processed are themselves first order in our accounting scheme.) While the 
effect of lensing on the CMB perturbation spectrum has been somewhat controversial, it 
can play an observable role on small angular scales. The fact that this second-order effect 
may be significant can be thought of as a consequence of the fact, noted in the Introduc- 
tion, that the existence of large distance scales in the problem can enhance higher-order 
effects; in this case the transverse deflection, given approximately by the product of the 
distance travelled times the lens angle, builds up as the photon travels along its trajectory. 

The term I'YY)k^^^'^di{(f) + t)^ is the longitudinal equivalent of the transverse lensing 
term. It arises from the time delay effect of the lenses, which alters the spacelike distance 
between the observer and the point where the photon path intersects the surface of emis- 
sion. The qualitative effect of this term is similar to that of the transverse lensing term, 
although its magnitude is expected to be smaller; for typical lens systems, the longitudinal 
deflection is smaller than the transverse deflection by a factor proportional to the lens 
angle (i.e., by several orders of magnitude). 

The term x^^^^do{(l) + i^)£ arises because the difference in affine parameter between 
observer and surface of last scattering differs along the true and background paths. It 
is similar in structure to the effects discussed in the previous two paragraphs, but is 
presumably smaller since the time derivatives of the potentials are typically smaller than 
the spatial derivatives. 
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Finally the integral hi^s) contains three terms. The first two terms appear small; 
they are integrals of products of two small quantities, and furthermore contain time deriva- 
tives which are typically suppressed with respect to spatial derivatives. The third term 
represents a correction to the ISW effect, taking into account that the perturbations along 
the first-order path differ from those along the background path. In cold dark matter mod- 
els with fl = 1 and adiabatic density perturbations, the ISW effect itself is smaller than the 
conventional Sachs- Wolfe term [8, 37-48], and the correction described here is presumably 
smaller still; nevertheless, it is possible that observations of the CMB will reach a level of 
precision at which this term should be taken into account. Moreover, time derivatives and 
the ISW effect can be important in models of structure formation based on topological 
defects [12-15], open universe models [49-54], and models with an appreciable cosmological 
constant [55, 56]. In these types of universes the new terms represented by l2{^s) could 
play a role analogous to that of gravitational lensing in adiabatic CDM models. 

V. Conclusions 

We have computed the anisotropy induced in the cosmic microwave background, due 
to gravitational effects, to second order in a given metric perturbation. For an arbitrary 
perturbation, our results are given by the basic equation (2.16) plus the solutions (3.13) 
for the perturbed geodesic and wavevector, where the forcing vectors to first and second 
order are given by (3.5). In the case of scalar perturbations in the longitudinal gauge, 
these results may be combined into the single compact formula (4.12). 

Our results are reassuring for studies to date of CMB anisotropy, in that they do not 
reveal any new effects which are likely to dominate the anisotropy spectrum on any scale. 
An informal examination of our final expressions indicates that the effect most likely to 
be observable is that due to (transverse) gravitational lensing, which has already been the 
subject of some attention in the literature. As both theoretical and observational studies 
of the CMB increase in accuracy and sophistication, however, we feel it is important to 
know the precise form of the effects we have explored. 

With the basic framework in hand, there is clearly room for future work along these 
lines. One direction would be to investigate a wider class of perturbations (i.e., vector and 
tensor modes), as well as to study carefully the second-order metric perturbation itself. An 
equally important task is to examine the effects we have described more quantitatively, in 
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the context of a specific and detailed model of structure formation; only then could we be 
completely confident in our understanding of the role played by second-order perturbations 
in CMB anisotropy. 
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Figure Caption 

Figure One. This figure shows the observer at location x'^(Ac)), the hypersurface of last 
scattering at r] = 'qg, and various paths connecting the two. The true geodesic in the per- 
turbed metric is x^{X), while the background geodesic is x^^^^{X). Adding the deviation 
vectors x^^^^{X) and x^'^^^{X) to the background path yields increasingly accurate approx- 
imations to the true path. The spacelike deviation vectors p*^"^* are to be distinguished 
from the a;^")'^, since the latter generally do not lie in hypersurfaces of constant conformal 
time. 
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